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Abstract 

We discuss azimuthal correlations in dAu collisions at different rapidities and centralities and argue that exper- 
imentally observed depletion of the back-to-back bump can be quantitatively explained by gluon saturation in the 
Color Glass Condensate of the Gold nucleus. 
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1. Introduction 

In a pioneering paper ||T| it was proposed to study the azimuthal correlations of hadrons produced at large rapidity 
separation Ay » 1 . The idea is that such correlations are mediated by the BFKL Pomeron. Therefore, unlike hadron 
production in hard collisions, where there is strong back-to-back correlation at opening azimuthal angle - n, 
correlations in the CGC should be significantly reduced. It has been suggested in ||2l that correlations at small Ay 
in the forward direction can also be used to study CGC. Indeed, forward direction corresponds to small x of nucleus 
where the CGC effects are strongest. They reduce both single and double inclusive hadron production and thus back- 
to-back correlations are suppressed. Unfortunately, there is a technical problem: the relevant scattering amplitudes 
are well-known in the so-called Multi-Regge-Kinematics Ay » 1, which is not applicable in this case. One therefore 
has to rely on phenomenological models, which offer descriptions that are analytically accurate only in parts of the 
interesting kinematic region. There are two such approaches: one that is based on the dipole model IS [31 and another 
one that is based on the fcj- -factorization |4|. 

The present calculation, based on 'A;7--factorization', assumes that 2 — > n process and the two-point correlation 
functions of CGC fields can be factored out. In this approximation, the 2 — > 4 amplitudes were calculated for an 
arbitrary Ay (quasi multi-Regge kinematics, QMRK) in |I5] |6l H ISl for gg — > ggqq and in ||8][T0l[lT] for gg gggg 
processes. Although generally A;?- -factorization fails in the gluon saturation region, there are valid reasons to believe 
that it provides a reasonable approximation of the observed quantities. Indeed, it was proved that A: j- -factorization 
provides the exact result for the cross section for single inclusive gluon production in the leading logarithmic approx- 
imation (LLA) Q [13.1 (though there is a subtlety in the definition of the unintegrated gluon distribution (f lil3lfT2l "). 
Although Atj -factorization fails for the double-inclusive heavy quark production, the deviation from the exact results 
is not large at RHIC energies 1 14|. At transverse momenta of produced particles much larger than Q,, A:^ -factorization 
rapidly converges to the exact results. There are also numerous indications that A:^ -factorization is phenomenologically 
reliable (see ||4| for examples). 
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2. Correlations at l3'r - 3'aI ^ 1 

First, we would like to consider correlations at small rapidity separations. Azimuthal correlation function is 
defined as 

C(A0)=^^, (1) 

where dN/d(A(p) is the number of pairs produced in the given opening angle and A^trig is the number of trigger 
particles. The number of pairs is given by 



dN 
d(A(p) 



2n [dkrkT [dyr [ dk^h [ dyJ '^^'^'^ ""^"^ i ] (2) 

J ^^J ^ J ^Xd^kTdjT d^kAdyA d^kjdyj d'^kAdjAj 



where kj- and yj are the transverse momentum and rapidity of the trigger particle and k^ and yA are the transverse 

momentum and rapidity of the associate one. We denote kj - -^k^ etc. throughout this paper. The first term on 
the r.h.s. of (j2| corresponds to gluon production in two different sub-collisions (i.e. at different impact parameters) 
and therefore gives a constant contribution to the correlation function, whereas the second term on the rh.s. describes 
production of two particles in the same sub-collision. The number of the trigger particles is given by 



dNtrig 

d^kjdyj 



A^trig = 2;r J dkjkj J dyr ■ (3) 



Expression for the single inclusive gluon cross section is well-known (see e.g. IfTSl ). The corresponding multiplicity 
reads 

dN 2 ffj 1 r T 9 9 

dyay^cTT^k^j' ^.^«(--^i)^^(-'(k-q,) )■ (4) 

In the center-of-mass frame x± = expliy) . Equation (4ji is derived in multi-Regge kinematics (MRK) x± <K 1. 
The correlated part of double-inclusive parton multiplicity is given by 

dNcon Nco] C d^qi f d^qi , , ^ 

5 (q, +%-\t- \a) 



r d'qi r 

d^krdyrd^kAdyA n^CfSj_ J q\ J q\ 

XipD(xx,q\)ipA(xi,ql)J[{c\^,(\^,'kj,\A^yT-yA), (5) 

where xi_2 - (kje'^^'^ + kAs'^^'^)/ yfs. The amplitude was computed in the quasi-multi-Regge-kinematics (QMRK) 
in ||8]|9][T0l and recently re-derived in fll] (the gg — » ggqq part was calculated before in ||5]|6l|7l|). In QMRK one 
assumes that xi,X2 <K 1, but Ay is finite. Explicit expression for can be found in ifTOl . 

For numerical calculations we need a model for the unintegrated gluon distribution function (p. In spirit of the 
KLN model l,15J we write 

= i 

where the saturation scale of nucleus is = ^'^^G^p. with Q^^^ the saturation scale of proton fixed by fits of the DIS 
data. The coupling constant is fixed at a, = 0.3. 

It has been pointed out in [ 10| that due to 1 — > 2 gluon splittings the double-inclusive cross section has a collinear 
singularity at 5 — > 0, i.e. it is proportional to [(Ay)^ -i- (A0)^]"'. Such singularities are usually cured at higher orders 
of perturbation theory. Additional contributions to the small angle correlations arise from various soft processes 
including resonance decays, hadronization, HBT correlations etc. Because the small angle correlations are beyond the 
focus of the present paper we simply regulate it by imposing a cutoff on the minimal possible value of the invariant 
mass s. This is done by redefining the amplitude as ^ ^ ^s/(n^ + s). For each kinematic region, parameter /u is 
fixed in such a way as to reproduce the value of the correlation function in pp collisions at zero opening angle A(p = 0. 

A;7-factorization is known to give results that are in qualitative agreement with a more accurate approaches, but 
miss the overall normalization. Therefore, in order to correct the overall normalization of the cross sections we 
multiply the single inclusive cross section Q by a constant Ki and the double-inclusive one (jSj) by a different constant 
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K2 lfT6l[T7ll . The correlation function C depends on both Ki and A'2. However, the difference Ca - C(A(p) - C(A0o) 
depends only on the ratio K2/K1. We choose A^o in such a way that C(A<po) is the minimum of the correlation 
function. This is analogous to the experimental procedure of removing the pedestal lilSi . The overall normalization 
of the correlation function K2/K1 - which is the only essential free parameter of our model - is fixed to reproduce the 
height of the correlation function in pp collisions. 

The results of the numerical calculations are shown in Fig. [TJ|3] In these figures we observe suppression of the 
bak-to-back correlation in dAu as compared to the bak-to-back correlation in pp, in agreement with the experimental 
data. In Fig. [3] we also see the depletion of the back-to-back correlation as a function of centrality. Note, that at the 
time of publication the precise centrality classes of the data shown in the lower row of Fig. [3] were not known. 




Figure 1: CoiTelation function at central rapidity. Kinematic region is 4 < pj- < 6, 2 < < (all momenta are in GeV), yr = 3.1, yA = 3. Left 
(right) panel: minbias pp (dAu) collisions. Data from |18| . 





Figure 2: Correlation function at forward rapidities. Kinematic region is pr > 4, 1.5 < Pa < Pt (all momenta are in GeV), yr = 3.1, va = 3. Left 
(right) panel: the minbias pp (dAu) collisions. Data from 1 19|. 



In addition to gg — > gggg and gg — » ggqq processes that we took into account in this section, production of valence 
quark of deuteron gqi, — > gq^gg gives a sizable contribution at forward rapidities due to moderate value of x associated 
with deuteron (x ^ 0.2 for pr = 2 GeV at y = 3). Contribution of this process to azimuthal correlations was analyzed 
in ||2l in the framework of the dipole model in MRK. However, the corresponding expression in Atj -factorization in 
QMRK is presently unknown thus preventing us from taking it into account in our calculation. In-spite of this we 
believe that the general structure of the correlation function as well as its centrality dependence is not strongly afiected 
by the valence quark contribution. We plan to address this problem elsewhere. 
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Figure 3: Correlation function at forward rapidities. Kinematic region is pr > 2, \.5 < Pa < Pt (all momenta are in GeV), yr = 3.1, Va = 3. 
Upper left (right) panel: minbias pp (dAu) collisions. Lower left (right) panel: peripheral (central) dAu collisions. Note: centrality of the theoretical 
calculation may not coincide with the centrality of the data (the former is not yet known at the time of publication). Data from 1 191 . 
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